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The problem of controlling dynamics of non- 
linear dynamical systems by stabilization of ex- 
isting unstable periodical orbits divides into two 
tasks. The first of them is picking out and 
analysing the periodical states. For this purpose 
the algorithm of Lathrop and Kostelich and algo- 
rithm of Schmelcher and Diakonos (SD method) 
are often used. The second task is stabilization 
of unstable orbit. In distributed systems the 
role of unstable orbits is performed by spatio- 
temporal periodical states. In our work we pro- 
pose methods of picking out unstable periodical 
spatio-temporal states with the help of modified 
Kostelich and SD methods for distributed beam- 
plasma systems. Stabilization is realized with 
the help of continuous delayed feedback, by the 
method based on Pyragas works. The most at- 
tractive feature of the proposed method is that 
the continuous control signal is given to one of the 
boundaries of the system what makes this method 
convenient to use in practice, for example, for mi- 
crowave beam-plasma systems. 



Introduction 

Practical use of unstable periodic orbits of chaotic at- 
tractor for chaos controlling has been discussed in a num- 
ber of papers since the the work of E. Ott, C. Grebogy, 
J.Yorke fil has been published. In the works 0, lEIiL 
IE IE IE HTE1 it has been shown that unstable orbits of 
chaotic attractor of a dynamical system with few degrees 
of freedom can be used to control complex behaviour of 
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the system, and the stabilization of the unstable periodic 
orbit needs small perturbation of the system dynamics. 
The work p| offers an algorithm of chaos control, defining 
the change of the system parameter in the dependence of 
the distance between the current state of the system and 
the wanted orbit. The value and sign of the parameter 
change can be discussed as the signal of some feedback. 
This method can also be applied to flow systems. 

As it has been shown in a number of works 0, fill IT3 . 
113j the OGY algorithm is insensitive to noise and inexact 
knowledge of the system state and it can be applied to 
the systems with finite, but large degrees of freedom [l4]| . 

But it seems rather hard to apply the OGY algorithm 
to control the dynamics of different distributed chaotic 
systems, including one we are discussing in this work, be- 
cause it demands the explicit knowledge of system state 
and the quick change of control parameter along the 
whole system space. It is nearly impossible to realize 
in real microwave beam-plasma systems working in the 
frequency rang 0.1 -j- 100 GHz. 

For such systems it seems more prospective to use the 
method of stabilization of the unstable periodic states 
suggested by K. Pyragas [l^. In this work the contin- 
uous |l6j . and not discrete quickly changing feedback is 
used, what makes this algorithm appropriate for the dis- 
tributed microwave systems. In the Pyragas scheme the 
system is synchronized with its own state taken one orbit 
period earlier, by continuous change of control parame- 
ter e(t) = 7(£(i) — £(t — T)), where is the analyzed 
variable, 7 is the feedback coefficient and T is the or- 
bit period. When the stabilization of the unstable orbit 
takes place, the feedback signal has the order of the noise 
level. It is important that all the information demanded 
for stabilization, except for T and 7, is contained in the 
time series of the system dynamics £(i), i.e. is defined au- 
tomatically in real time. Applications of such controlling 
scheme were theoretically and practically investigated in 
different systems and models @, Hfl UJl HE liil including 
laser physics [2E El E3 1 models of geophysical processes 
[2^ and the processes of "reaction-diffusion" j24[ . 
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We must mark especially that this scheme with con- 
tinuous feedback has been used for controlling spatio- 
temporal chaos in distributed systems (see for exam- 
ple works pil I2H I2IH devoted to stabilization of two- 
dimensional structures in the 2D distributed chaotic sys- 
tem, described by partial derivative equations; works 
[27l I2H |29| . in which chaos control in the lattices of 
coupled maps has been considered; the works analyzing 
such sample equations as the complex Ginzburg-Landau 
equation in non-linear dynamics |3Ct ]3l| or the Swift- 
Hohenberg equation of laser physics [3l 1331 '). Among 
such investigations we must mark wor ks | 24l l3lT | espe- 
cially. Thus, in the work of Boccaletti [3l| the problem 
of chaos control in the Ginzburg-Landau model with the 
use of a large (but finite) number of small local perturba- 
tions has been investigated. This method allows to sta- 
bilize the unstable structures of chaotic spatio-temporal 
dynamics and to synchronize two chaotic states. Frances- 
chini and the others in work 24] studied the problem of 
controlling chaotic generation of impulses in the model 
of uniformly coupled "reaction-diffusion" systems used 
for describing the processes of charge transfer in bistable 
semiconductor. The controlling feedback signal in math- 
ematical model offered in work [24| is supposed to be 
spatially homogeneous, i.e. the perturbation carried into 
every point of the system space is the same. 

Some works describing chaos control in microwave 
beam-plasma systems must also be mentioned. The sta- 
bilization of unstable periodic orbits in a fluid model of 
Pierce diode with the help of the OGY algorithm has 
been discussed in work [34| . Work [35| discusses the prob- 
lem of suppression of the regime of the virtual cathode 
in Pierce diode |36l 1371 ] with the help of delayed feedback 
using a numerical method of large particles. In our work 
[3l| the influence of the external feedback on the chaotic 
oscillations in Pierce diode has been examined. 

It is important to note, that in the majority of works 
devoted to chaos control in the distributed systems, the 
controlling feedback signal effects at every point of the 
system space. Such distributed spatio-temporal control 
can be used, for example, in optical systems [39ll40l l4l | . 
But application of such methods to beam-plasma sys- 
tems causes several difficulties. In the present work we 
propose a method of chaos control based on the continu- 
ous feedback scheme, where the feedback signal is given 
to one of the boundaries of the distributed system. As 
a sample this method is applied to widely known simple 
distributed beam-plasma system (fluid model of Pierce 
diode) 

As it has already been mentioned above, the knowledge 
of the positions of the unstable orbits in the phase space 
and their periods is quite important for controlling chaos 
in the discrete and flow systems with a small number of 
degrees of freedom. We can say that in distributed sys- 
tems the unstable periodic spatio-temporal states play 
the roles of unstable orbits. Therefore, for the efficient 
chaos controlling we need to find such unstable periodic 
states. Accordingly, we propose a simple method of de- 



tecting of unstable periodic spatio-temporal states and 
their further stabilization. 

The structure of the paper is following: In Section[l]the 
fluid model of Pierce diode is discussed. In SectionlTTIsta- 
bilization of the unstable homogeneous equilibrium state 
is described. In Section II I II the spatio-temporal data 
series of electron beam dynamics are analyzed and the 
unstable periodic spatio-temporal states of the chaotic 
dynamics are carried out. And, in Section Hvl we discuss 
the stabilization of the unstable periodic spatio-temporal 
states (unstable states) derived in Section ITTT1 with the 
help of feedback of different type. In conclusion we sum- 
marize the main results discussed in our paper. 



I. GENERAL FORMALISM 

Pierce diode [3(| |3?J is one of the simplest beam- 
plasma systems demonstrating complex chaotic dynamics 
IH S El El 0- It consists of two plane parallel infi- 
nite grids pierced by the monoenergetic (at the entrance) 
electron beam (see Fig.^J. The grids are grounded and 
the distance between them is L. The entrance space 
charge density po and velocity vq are maintained con- 
stant. The space between the grids is evenly filled by the 
neutralizing ions with density \pi/po\ = 1. The dynam- 
ics of this system is defined by the only parameter, the 
so-called Pierce parameter 



uj p L/vo, 



(1) 



where u p is the plasma frequency of the electron beam. 
With a > ir in the system, the so-called Pierce instability 
[3^. l37j develops, which leads to the appearance of the 
virtual cathode. At the same time, with a ~ 3tt, the 
instability is limited by non-linearity and the regime of 
complete passing of the electron beam through the diode 
space can be observed. In this case the system can be 
described by the fluid equations: 
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with the boundary conditions: 

«(0,*) = 1, p(0,f) = l, <p(p,t)=<p(l,t) = 0. (5) 

In equations the non-dimensional variables 

(space charge potential <p, density p, velocity v, space 
coordinate x and time t) are used. They are related to 
the corresponding dimensional variables as follows: 



^ = (vl/v)^, E> = (v 2 /L V )E, 
p' = pop, v' — vqv, x' — Lx, t' — (L/vo)t, 



(6) 
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FIG. 1: Schematic diagram of Pierce diode with a feedback 
line 

where the dotted symbols correspond to the dimensional 
values, rj is the specific electron charge, vq and po are the 
non-perturbed velocity and density of the electron beam, 
L is the length of the diode space. 

Equations © and © are numerically integrated with 
the help of one-step explicit two-level scheme with up- 
stream differences [U EH m and Poisson equation Q 
is solved by the method of the error vector propagation 
|45|. The time and space integration steps have been 
taken as Ax = 0.0005, and At = 0.0001, respectively. 

In works |37l l42l |43 | it was shown, that in the nar- 
row range a ~ 37r without the virtual cathode, one can 
observe chaotic oscillations in the system. With a de- 
crease of a from 2.887T to 2.867T the model demonstrates 
a transition from periodic oscillations to chaos via a pe- 
riod doubling cascade. Just below the critical value of a 
the system demonstrates weak chaotic oscillations with 
the explicit time scale. Below we will call this regime 
"bond chaos" because the attractor in the reconstructed 
phase space looks like a narrow bond on which the phase 
paths are situated. With the further decrease of a the ex- 
plicit time scale vanishes and the phase picture reminds 
of spiral twisting from one point. The power spectrum 
of the system oscillations is more complex in comparison 
with the first case. This type of chaotic behaviour we 
call "spiral chaos" . 

The stationary homogeneous equilibrium state of the 
electron beam in Pierce diode is characterized by the fol- 
lowing distribution of the system variables 

v(x) = 1.0, p(x) = 1.0, <p(x) = 0.0, (7) 

that can be found from solutions of stationary equations 
derived from non-stationary fluid equations (J2J Q in as- 
sumption that d/dt = 0. 

Homogeneous equilibrium state Q is stable with a < 
7r, and loses its stability at a > ir as it has been shown by 
J.Pierce [3(j. In work [3^| the significance of the unsta- 
ble a decrease of a or with the influence of the external 
delayed feedback is shown. 



II. STABILIZATION OF UNSTABLE 
EQUILIBRIUM SPATIAL STATE WITH THE 
HELP OF THE CONTINUOUS DELAYED 
FEEDBACK 

Follo wing the works mentioned above 0, IrH fTsL Il9l 
H3, El El El devoted to chaos control in finite- 
dimensional systems, we now describe the scheme of the 
continuous delayed feedback for stabilization of the un- 
stable equilibrium state in Pierce diode. 

The continuous delayed feedback is realized by chang- 
ing the potential of the right boundary of the system: 

tp(x = 1.0,*) = /fb(*) = K(p(x Rx ,t)-p(x Rx ,t-d)), (8) 

where K is the feedback coefficient, d is the delay time. 
Variable p(xa x ,t) in (jSJ is the space charge density in a 
fixed point x = xa x of diode space (in the work we use 
Xfix = 0.2). When the stabilization regime takes place 
and the system is exactly in the equilibrium state, the 
feedback signal /fb could by compared with the noise 
level. 

Practically, this scheme of delaying feedback can be re- 
alized, for example, with the help of delay lines on mag- 
netostatic waves 47, 481 or acoustic waves [Hill. Such 
method would allow to select the required delay time. 
The signal can be picked out from the system using a 
probe. 

so, unlike works mElElElElEllIlliiEi: 

where controlling the chaotic dynamics of the distributed 
systems needs the effect of the feedback signal upon 
the whole interaction interval, in our scheme the signal 
of the continuous delayed feedback determines only the 
change of the boundary conditions. It is more convenient 
to realize in experimental microwave devices. Then, in 
real microwave systems working in the frequency range 
0.1 -i- 100 GHz, it is nearly impossible to realize the OGY 
algorithm because it demands very quick changing of con- 
trolling parameter. 

It must be marked that we use the term continuous 
feedback in analogy with works [l5l Il6| . This term 
emphasizes that in our scheme the controlling signal is 
changed continuously in contrast to the schemes based 
on the OGY algorithm pj. It is obvious that the control- 
ling scheme (JSJ discussed in our work is more convenient 
in practical use for controlling the dynamics of the dis- 
tributed microwave systems. 

For purpose of stabilization of the unstable equilibrium 
state the delay time d of the feedback line © must be 
small enough: d < T p /2, where T p — tt/lo p is the typical 
time scale in the autonomous distributed system |37| . 

Results of stabilization of the system ©-(0) with the 
help of the continuous feedback © are presented in 
Fig.El The autonomous system displays chaotic oscilla- 
tions with a large amplitude. But after switching on the 
continuous feedback one can observe a fast decrease of the 
amplitude of oscillations which leads to the stabilization 
of the unstable state Q . After a short transient process 
the controlling signal in the feedback line becomes rather 
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FIG. 2: Time series p(a:fi x , t) in autonomous system (top pic- 
ture), oscillations pa(xa x ,t) in stabilized system (middle pic- 
ture) and the signal in feedback line fib{t) (bottom picture) 
for the spiral chaos regime (a = 2.8577r). The feedback pa- 
rameters are K — 0.8, d — 0.15. The moment of appearance 
of the delayed feedback is marked by an arrow and dotted 
line. 



small in comparison with the signal before stabilization 
(near 0.01%). This means that the regime of chaos con- 
trol with the help of small controlling signal of continuous 
feedback takes place in the system. 

It is important to find the regions of feedback param- 
eters K and d where chaos controlling is possible. Fig-El 
shows areas of the bond chaos regime (a = 2.8627r) and 
spiral chaos regime (a = 2.85 77r). One can see that the 
regime maps {K 1 d) for both cases have qualitative sim- 
ilarities. For small K, the system displays chaotic oscil- 
lations similar to autonomous oscillations corresponding 
to this a value (region C in Fig.[3J|. With an increase 
of K one can observe the destruction of chaotic oscilla- 
tions and appearance of periodical ones (region P). In 
this case the feedback signal is not small (its amplitude 
has the same order of value as /fb before switching on 
the feedback) hence we can not regard this regime as 
the controlling chaos one. With further growth of K the 
controlling chaos regime takes place (region S) . The un- 
stable homogeneous state is stabilized, and the system 
dynamics in this case is shown in Fig.[3 

The width of region of stabilization S depends notice- 
ably on the delay time d. The threshold values d\, di, 
determine the range d £ (di,e?2) where controlling pa- 
rameter values allowing to obtain controlling chaos in 
the distributed beam-plasma system are possible. 

With the further increase of d and K, reflection and 
overtaking of electrons appear in the electron beam (re- 
gion V in Fig.[3J), and therefore the fluid equations ©-© 




FIG. 3: Regime maps of the stabilized system on the plane 
of parameters (K, d) for cases of (a) bond chaos (a = 2.8627r) 
and (b) spiral chaos (a = 2.8577r). The region of stabilization 
of the unstable equilibrium state is hatched 



become unfit (see |5l|,|52j for detail). 

It is useful to note that in spiral chaos regime the region 
S in the parameter map is narrower than that in bond 
chaos regime (compare Fig.^Ji and [210. 

We can analyse the stability of the discussed equilib- 
rium state by calculating the maximum Lyapunov expo- 
nent A [l5| ■ For the feedback coefficient K — the values 
of Lyapunov exponents of our system are equal to those 
of the autonomous system. 

Fig.0] shows the dependence between the maximum 
Lyapunov exponent A and the feedback coefficient K for 
two different chaotic regimes. This dependencies deter- 
mine the boundary of adaptability of our control method. 
The suppression of the chaotic dynamics is possible only 
for the values of K corresponding to A(K) < 0. These 
are regions P and S in Fig.|3J The value of K for which 
stabilization of unstable state takes place (region S in 
Fig-EJ) is marked by an arrow in Fig.0] One can see that 
the stabilization regime is preceded by the increase of 
Lyapunov exponent, and beyond the boundary of stabi- 
lization K = K c the value A(K) decreases linearly. Such 
behaviour is typical to both types of chaotic dynamics. 
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FIG. 4: Dependence of the maximum Lyapunov exponent A 
upon the feedback coefficient K. The dotted line corresponds 
to the bond chaos regime (a = 2.8627r), the solid line — to the 
spiral chaos regime (a = 2.8577r). The arrow marks the values 
of K c , corresponding to the onset of stabilization, d — 0.14 



III. PICKING OUT UNSTABLE PERIODICAL 
SPATIO-TEMPORAL STATES OF SPATIALLY 
EXTENDED SYSTEM DYNAMICS 

In this section we discuss the unstable periodic spatio- 
temporal states of distributed active medium, analogous 
to the unstable periodic orbits of a few-dimensional sys- 
tems, which play an important role in the complex be- 
havior of non-linear systems l53|. in particular, in the 
case of chaotic synchronization [5J, |5j| llil |53, IH, |5!j • 

The initial information about the set of unstable pe- 
riodic orbits can be obtained from the histograms de- 
scribing the frequency of system returning to the vicinity 
of the orbits. This method was offered by D.P. Lathrop 
E.J. Kostelich We choose a phase point Ri on the 
attractor. If this point is close to an unstable cycle with 
a period T, then the phase trajectory passing the points 
Rj+i, Rj+2 7 • ■ ■ i R-fc ; will come near the initial state with 
some precision s > 0: 



R? — R 



i+m I 



(9) 



where m = T / At is the orbit period in discrete time 
unit. Then the distribution of return time (histogram) is 
made. By this histogram it is easy to carry out the typical 
periods of corresponding unstable states, and then to find 
the unstable cycles. 

For the study of unstable spatio-temporal states of the 
dynamics of the explored distributed system we inves- 
tigate the time oscillations of the space charge density 
p(xo,t), obtained from the fixed points of the interac- 
tion space xq. Then, from the time series p(xo,t), using 
Takens' method [6(J the attractors of system dynamics 
are restored in pseudo-phase space R n . Using Lathrop 
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FIG. 5: Histograms showing the number of visits of phase tra- 
jectory to the fixed attractor points, associated with different 
recurrence time: (a) a — 2.86l7r; (6) a — 2.8577T 



& Kostelich algorithm we pick out the periodical orbits 
in R n and then hence find the spatio-temporal periodical 
states of the distributed system. 

In Fig. the histograms of return time obtained for 
different values of Pierce parameter are presented. One 
can see that in the bond chaos regime (Fig. [SJ one of the 
orbits dominate upon the others i.e. it is more frequently 
visited by phase point. 

For spiral chaos regime (Fig. 0>) the set of unstable 
orbits is more complex, and the different unstable peri- 
odical states are visited more evenly and frequently. 

Histograms of recurrence time obtained for time series 
p(x,t), taken from different points of diode space x = Xg 
are very much alike. This is illustrated in Fig. It 
follows from the picture that the set of unstable spatio- 
temporal states is identical along the diode interval. A 
slight difference can be caused by inaccuracy in calculat- 
ing the recurrence time (see, for example, the unstable 
spatio-temporal states with T = 4.173 in Fig. |SJa and 
Fig- lit (marked by an arrow)). The latter allows to carry 
out and to analyse some characteristics of the unstable 
spatio-temporal states of the distributed system using 
scalar time series obtained from one fixed point Xfi X of 
interaction interval. 

In Fig.[7| the unstable periodical states of the dis- 
tributed system are shown. The three-dimensional pic- 



6 



N(T) 

120 - 



T=4.173 



80 



4G 



/ 



a 



i-i.jj|J u r .L A, ..I.I.J,.. t 



N(T) 

120 - 
80 
40 - 



20 40 60 80 100 





N(T) 

120 



80 



/ 



ImLL 



j^ljiaJ Urn 



T 



20 40 60 80 100 



40 








20 



-10 



60 



80 



100 



FIG. 6: Histograms of recurrence time for a — 2.857-7T ob- 
tained from the data taken from different points of diode 
space: (a) xo = 0.2, (b) xo = 0.5, (c) x — 0.8 



tures show the spatio-temporal oscillations of the space 
charge density p(x, t) which correspond to the unstable 
periodical orbit states with different period. The pictures 
are obtained for a = 2.8577T (spiral chaos). 

We must pay attention that the similar results has 
been derived when the unstable periodic spatio-temporal 
states by the method of P. Schmelcher F. Diakonos (SD- 
method) [jj 0, adapted for the analysis of distributed 
system. The algorithm is the follows. 

As in the case when the unstable periodic states 




FIG. 7: Unstable periodical orbits of different periods T in 
the regime of spiral chaos: (a) T = 4.173 and (6) T = 18.987 



are reconstructed from histograms, at first we deter- 
mine the system state vector R in phase space. As 
a state vector in this case a vector is taken com- 
bined from the density of charge taken from the 
different points of the interaction space: R(i) = 
{p(x = 0.25, t),p(x = 0.5, t),p(x = 0.75, t)} T . 

In reconstructed phase space a plane p(x = 0.25, t) — 
1.0 is taken, which is considered as Poincare section. Let 
us indicate as R„ the system state corresponding to the 
nth crossing between phase part and the selected section, 
(this crossing takes place at moment of time t n ). Then 
the description of the system can be made with the help 
of discrete map like 



Rn+i — G(R„), 



(10) 



where G(-) is the evolution operator. Obviously, it is im- 
possible to find the analytical form for the operator G, 
but numerical integration of the initial system of fluid 
equations can give us a consecution of values {R}„, de- 
rived from the map l|10f) . 

SD-method for picking out unstable periodic orbits 
supposes examination of modified model of a distributed 
system, described by following map 0: 



R 



n+l 



Rn 



AC [G(R„) - R„] . 



(11) 



In the works 0, 13 ^ is exactly shown that such modi- 
fied system in the case of analysis of discrete systems and 
the systems with few degrees of freedom allows to stabi- 
lize effectively the unstable periodical orbits of the initial 
system, which in the case of modified system transform 
from saddle into stable ones l|llf> . Application of the SD- 
method to the analysis of distributed system is also ef- 
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fective for picking out unstable periodic spatio-temporal 
states. 

In the equation A = 0.1 is the method constant 
and C is a matrix, which, following to work Qj was taken 
as 

c-(J°). (i*> 

Modified map Hll|) allows only to define unstable pe- 
riodic state of lowest period T\, which agrees with the 
single crossing of the selected Poincare section during 
time T\. To find the unstable states with higher period 
T p , when the phase path crosses the Poincare section for 
p times during the period T p , it is necessary to take a 
modified map like 

R-n+i = R n + AC [G P (R„) — R„] , (13) 

where G p (-) is p times iterated map G(-) (i.e. when nu- 
merical solution of the system of fluid equations is found 
and then the map is reconstructed 1|13[) it is necessary 
to take into account only the every pth crossing of the 
Poincare section by phase path, where p is obviously in- 
teger. 

So, by numerical integration of the map (f 1 31) with 
different values of p one can find the unstable spatio- 
temporal states that appear to be stable in the modified 
system, described by the map l|13|) . So, the result of nu- 
merical integration of modified system (fluid equations 
taking into account the procedure is a periodic so- 

lution, which agrees with the unstable periodic state of 
the initial system. 

The numerical integration of the modified system 
causes only one difficulty. When we are searching for 
the state R n +i at the moment t n+1 , then, according to 
formula H13f) we know only the coordinates of this state 
in Poincare section but we don't know the correspond- 
ing distribution of space charge density p(x,t n+1 ), ve- 
locity v(x,t n+1 ) of the electron beam and the potential 
4>{x,t n+1 ). To derive the mentioned space functions we 
use the following procedure. The system of fluid equa- 
tions describing Pierce diode is integrated till some state 
R s of the system would not be close to the selected state 
R„+i with some demanded precision: ||R, i+ i — R s || < 5, 
where 8 is taken as <5 = 10~ 3 . When this condition is 
fulfilled, the space functions corresponding to the state 
R s are considered as the space distributions R„+i and 
then the next iteration according to l|13|) takes place. 

As numerical modelling shows, the modified SD- 
method applied to distributed system is convergent and 
allows to find the demanded periodical time-space states. 
The convergence of the procedure is illustrated by Fig- EI 
which shows the depcndance of the space charge density 
p(x = 0.25) in the moment when the point in pseudo- 
phase space crosses the Poincare section upon the num- 
ber of iteration n of the SD-method when the unsta- 
ble state of period 1 is defined (p = 1). One can see 
clearly that the iteration process of SD-method con- 
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FIG. 8: Dependance of the space charge density p taken at 
the moment when the phase point passes the Poincare section, 
upon the number of iteration of SD-method when the orbit of 
period 1 (T = 4.2) is picked out. Pierce parameter a — 2.8587T 
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FIG. 9: Time series of the space charge density p(t) in the 
fixed point of the phase space x — 0.2, corresponding to the 
different unstable periodic spatio-temporal states picked out 
by SD-method with a = 2.858tt: 1 - T = 4.2 (p = 1), 2 - 
T = 9.1 (p = 2), 3 - T = 16.9 (p = 3), 1 - T = 18.9 (p = 4) 



verges to the value corresponding to the unstable time- 
periodical spatio-temporal state of the system. Fig. [5] 
shows the time series of p(t,x — 0.2), corresponding to 
unstable periodic states of different order p and period T 
with a — 2.8587T, picked by SD-method. We mark that 
the derived spatio-temporal distributions agrees very well 
with the distributions of unstable states derived earlier 
by analysing of the return time histograms. 

As quantitative characteristics of the derived unsta- 
ble orbits it is useful to analyze the maximum Lyapunov 
exponent A of every orbit. Calculation of this value is 
important for further stabilization of periodical spatio- 
temporal states. 

Values of maximum Lyapunov exponent A T calculated 
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Table. 1: Values of maximum Lyapunov exponent A T for 
unstable states with period T in the regime of spiral chaos 
(a = 2.857tt) 

Period of unstable Maximum Lyapunov 



state T 


exponent A T 


4.347 


0.854 


8.289 


1.194 


10.581 


0.698 


12.501 


0.629 


14.613 


0.337 


18.987 


0.266 


23.115 


0.186 


39.636 


0.094 


43.953 


0.083 


60.438 


0.053 



with the help of Benettin's algorithm [§lj for the most 
frequently visited orbits are presented in Table 1. 



IV. STABILIZATION OF THE UNSTABLE 
PERIODICAL STATES 
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As simplest scheme of stabilization of unstable period- 
ical orbits described in the previous section, we can take 
the scheme (JSJ where the feedback signal is formed as: 



tp(x = 1.0, t) 



/fb 



(*) 



FIG. 10: Spatio-temporal dynamics of Pierce diode in chaotic 
regime (a) and in the regime of stabilization of the periodic 
orbit T = 4.173 (b). The arrow and dotted line in b shows 
the moment of switching on the delayed feedback 



= K(p{ Xfix , t) - p{ Xfix , t - T k )) = K£(t). (14) 

Here K is again the feedback coefficient and 7\. is the 
delay time equal to the period of the fc-th unstable orbit. 
As in the case of stabilization of unstable equilibrium 
state, we choose the fixed point of the diode space £g x = 
0.2, from which the feedback signal is taken. 

The numerical modelling shows that this scheme is 
rather effective for stabilization of the unstable periodical 
state with the lowest period T\. Fig.lTfjIshows the spatio- 
temporal dynamics of the system in cases of autonomous 
oscillations and in the regime of stabilization (the value 
of the space charge density p(x, t) is shown by colour nu- 
ances scaling). One can easily observe that there are pe- 
culiarities of the complex spatio-temporal behaviour and 
transitions between chaotic and periodic regimes. It is 
obvious from the picture that the periodic dynamics anal- 
ogous to the system behaviour near the unstable state is 
stabilized during approximately 2-^3 time periods T\ (see 
Pig-Eh). 

We also study the peculiarities of the influence of de- 
layed feedback (|14f> upon the investigated system by ex- 
amining how the feedback amplitude K acts upon the 
maximum Lyapunov exponent A and the average of the 
feedback signal 

(£> = - f fl(t)dt, (15) 
Jo 



where r 3> d. The corresponding plots derived for d = T\ 
are presented in Fig. ^] 

When K is small the system demonstrates chaotic os- 
cillations similar to the oscillations in the autonomous 
regime (in Fig. 112b . one can see the time series of the 
space charge oscillations for x = 0.2 without feedback). 
With an increase of K the complexity of the oscilla- 
tions diminishes (the value of A decreases) and at the 
same time the amplitude of the feedback signal becomes 
smaller. In some region of K the feedback signal is close 
to zero and the maximum Lyapunov exponent A < (the 
grey area in Fig. lllfl . In this case the spatio-temporal dy- 
namics is the same as that of the unstable periodic state. 
Thus, it is the regime of chaos stabilization, illustrated 
by Fig. 112b. c. which shows the oscillations of the space 
charge density in the system and the feedback signal, 
correspondingly. 

Examining the graph A(K) we also must note the re- 
gions where the feedback signal (£) is not small but the 
value A < 0. These regimes are marked in Fig. ^2 as P 
and correspond to periodic oscillations near the unsta- 
ble orbit. The space charge oscillations and the feedback 
signal for this case are presented in Fig.ll2ti. e. We also 
pay attention to the fact, that with large K the system 
demonstrates the period doubling based on the regime 
P (region marked as P 2 in Fig. lll|) . For K > 0.15 in 
the system the growth of the oscillation amplitude takes 
place, and reflection of electrons appears. Thus, in this 
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FIG. 11: Dependence of the maximum Lyapunov exponent 
(solid line) and the mean feedback signal (doted line) upon 
the feedback coefficient K in the spiral chaos regime (a = 
2.8577t). Region C (grey colour) corresponds to the regime 
of stabilization of the unstable periodic state of period 1 T = 
4.173 (the feedback signal is close to null (£) « 0). Region P 
corresponds to periodic oscillations, when the feedback signal 
is large and the maximum Lyapunov exponent A < 0; region 
P 2 corresponds to the regime of period doubling 



region we can not use any more the fluid equations for 

i i I i 

the description of the system |5JJ, |52( . 

But the computer experiment shows that stabilization 
of the unstable orbits with higher periods > T\ = 
4.173 using feedback scheme is impossible. The anal- 
ysis of the chaos controlling schemes like (|14|) applied to 
the systems with few degrees of freedom, has shown j^l 
that the effective stabilization of the unstable periodic 
orbits is possible only when the values of the maximum 
Lyapunov exponent A and the orbit period r fulfill the 
following condition 



At < C, 



(16) 



where C is constant depending on the system. 

Thus, we must modify the controlling scheme l|14f> so 
that it would be possible to stabilize the orbits with 
higher period The scheme, which we consider, is 
a modification of the method worked out by K.Pyragas 
[To| . As it follows from the result of the numerical ex- 
periment [6^, in the few-dimensional non- linear systems 
this scheme allows to stabilize the unstable orbits for 
which the condition Ijlbl does not fulfill. The main idea 
of the method ||33| is the following: the feedback sig- 
nal is formed so that not only the system state at the 
moment (t — T&) influences upon it, as it was in the pre- 
vious scheme, but also the states at the moments of time 
(t — mTk) do, with some weight coefficients. Following 
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FIG. 12: Time series p(x^,i) of the autonomous system (a), 
time series p a {x^,i) in the stabilized system (K = 0.07, d = 
T — 4.173) (6), feedback signal fo,(t) (c); time series p p (xf LVi , t) 
and /fb(i) (d,e) in periodic oscillation regime (A < 0) (K = 
0.035, d = T = 4.173), a = 2.857tt (spiral chaos regime). The 
moment of switching on the feedback is marked by an arrow 
and a dotted line. 



work |2J] , we assume that the described scheme could be 
effective in the controlling dynamics of the distributed 
chaotic system. 

The continuous feedback for stabilization of the un- 
stable periodic states of higher period Tk > T\ can be 
described by the following expression 



ip{x = 1.0, t) = K£(t) 



= K^(l-r)J2 r m - X Mara*, t) ~ p(x &x , t - md))j , 

where the delay time d is taken equal to the period of the 
stabilized periodic orbit T^, M is large enough (M ^> 1), 
and r (0 < r < 1) characterizes the contribution of the 
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previous system's states into the feedback signal: small 
values of r correspond to the small significance of the 
previous states in the signal ifTTjl . and analogously, cor- 
respond to the large ones. The case r — corresponds 
to the simplest scheme (I14|) of continuous feedback dis- 
cussed above. Such feedback scheme Ijl7|l can be realized 
in practice, for example, by using a set of acoustic time 
delay lines each with its own delay time and transfer con- 
stant. As it follows from works [24], I6J, |65j , a possibility 
of stabilization of the orbits of highest period T k ljTo|) 
appears with the growth of r, for which the use of the 
scheme l|14[l is impossible. For example, in the work [fj^] 
it is shown that with the help of the scheme similar to 
(fTTjl it is possible to stabilize the unstable periodic states 
for which the following estimation condition is fulfilled: 

1 + r 

AT<C^—. (18) 
1 — r 

Using the scheme (|17|l for the stabilization of the peri- 
odic orbits which fulfills the condition (|18fl . the explicit 
knowledge of the orbit period Tk is needed. If the de- 
lay time d differs even slightly from the real period Tk, 
it is impossible to stabilize the unstable periodic state 
of the system, and in the feedback line one can observe 
periodic oscillations with the amplitude exceeding 
significantly the noise level in the system. As numerical 
analysis shows, stabilization is possible only when the 
orbit period is known with the precision e < 0.01%. 

For more accurate definition of the orbit periods Tk 
the following method presented in work |66( for the sys- 
tems with few degrees of freedom and approbated for 
distributed systems in the work [24| is used. 

If the delay time d is not explicitly "tuned" to the orbit 
period Tk then the maximum Lyapunov exponent is neg- 
ative and in the feedback line one can observe the oscilla- 
tions with the amplitude exceeding significantly the noise 
level of the system and the "base" period 8^^, which 
depends upon K and d (see. Fig.ll2li.e). In work 
the analytical dependence between the unknown explicit 
period Tk of the unstable state, feedback parameters K, 
d and the period G has been found as follows: 

B(K,d) =T k + — —(d-Tk) + 0((d-T) 2 ). (19) 
K-g 

Here g is the unknown parameter which is defined by the 
type of the non-linear dynamical system and depends on 
the form of the delayed feedback. Thus, in expression 
(|19|l there are two unknown values: the period Tk of the 
unstable periodic state and the parameter g. 

Let us take two sets of values {Kx,d\), (K 2 , g^), where 
the values of dip are chosen close to the period of the 
stabilized orbit, and then find the corresponding Gi and 
02- Then we find numerically the solution of the sys- 
tem of two non-linear equations l|19|) and determine the 



values g and Tk- Then we repeat the procedure taking 
the newly defined Tk as d. After a number of such it- 
erations the period of the unstable orbit Tk, would be 
defined with absolutely satisfactory precision and then 
the scheme (|17|) could be used for stabilization of this 
state. 

Numerical analysis shows that with the help of this 
method it is possible to stabilize the orbits with the pe- 
riod T < 25.0. The unstable states with higher period 
turned out to be impossible to stabilize, though in the 
system the periodic oscillations can be observed, but the 
form of these oscillations is not close to the concerned un- 
stable state and the feedback signal is not small. With 
the growth of the period T of the stabilized orbit it was 
necessary to enlarge the parameters M and r to achieve 
the effect of controlling chaos. Simultaneously, the range 
of value of the feedback coefficient K, for which stabi- 
lization of the unstable state is possible, narrows. 

Stabilization of the orbits of highest period is illus- 
trated by Fig. The feedback scheme parameters M, 
r and K arc defined in the outline. The explicit values 
of the periods T are defined according to the method de- 
scribed above. In Fig. ^|the orbits of the period T f» 8, 
T » 16 and T w 23 are shown. The latter is the unstable 
orbit with the highest period which appears to be possi- 
ble to stabilize in spiral chaos regime (a — 2.85 77r) with 
the help of the discussed scheme. 



Conclusion 

In this work the method of controlling complex chaotic 
dynamics of the spatially distributed active medium 
"electron beam with overcritical current in Pierce diode" 
is discussed. The described method is based on the ideas 
of controlling chaos in non-linear systems with few de- 
grees of freedom. The scheme of continuous delayed feed- 
back, which is used for controlling, allows to stabilize the 
unstable equilibrium state of the distributed system and 
the unstable periodic spatio-temporal states analogous 
to the unstable periodic orbits of the chaotic attractor in 
the systems with few degrees of freedom. 
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FIG. 13: Stabilization of the unstable states of high period: (a) T = 8.3829 (M = 3, K = 0.03, r = 0.2); (b) T = 16.6582 
(M = 12, K = 0.012, r = 0.38); (c) T = 23.1286 (M = 20, K = 0.008, r = 0.56). The pictures (from left to right) present: 
the time series of feedback signal /ft,, the power spectrum of the oscillations in stabilized system and the phase portrait for the 
space charge density oscillations p(xa^ = 0.2, t). The feedback is switched on at the moment of time t = 50.0. Delay time of 
time-delayed embedding 5 = 1 
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